Fidelities for transformations of unknown quantum states 
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We present a general theoretical formalism to compute the fidelity of transformations of unknown 
quantum states, and we apply our theory to Gaussian transformations of continuous variable quan- 
tum systems. For the case of a Gaussian distribution of displaced coherent states, the theory is 
readily tractable by a covariance matrix formalism, and a wider class of states, exemplified by Fock 
states, can be treated efnciently by the Wigner function formalism. Given the distribution of input 
states, the optimum feed back gain is identified, and analytical results for the fidelities are presented 
for recently implemented teleportation and memory storage protocols for continuous variables, 
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I. INTRODUCTION 

In a genèric scenario for quantum state transforma- 
tions, a protocol is applied to an unknown input quan- 
tum state taken from a family of states with a certain 
probability distribution. The quality of the protocol is 
quantified by a fidelity measure, which in a natural sense 
extracts the average overlap between the state obtained 
through use of the protocol and the state expected under 
ideal circumstances. 

In laboratory experiments, one may delegate the han- 
dling of the initial random state preparation and the fi- 
nal examination of the output to an independent person 
or device, Victor, (and appeal that Victor is not leak- 
ing information to the other experimenters). Theoretical 
physics works differently in the sense that what is spec- 
ified at one point in the theory is specified throughout, 
and quantum physics adds the further aspect that the 
quantum state is a state describing our knowledge about 
the systcm, influcnccd by any knowledge that comes to 
our mind be it in the form of measurement outcomes or 
information about the preparation procedure for the sys- 
tem - we cannot know something about a system and at 
the same time describe it by a state vector or density 
operator that is independent of this knowledge. 

It is the purpose of this paper to present a practi- 
cally useful theory to determine fidelities without leaving 
doubts about the correct handling of what is known and 
what is not known about the input states. In Sec. II, we 
introduce a theoretical formalism which represents Vic- 
tor both at the stage of preparation of an unknown input 
state and at the examination of the output. In Sec. III, 
we consider the case of continuous variable quantum sys- 
tems, and we show that a hybrid quantum and classical 
Wigner function can be applied in an explicit calculation 
of the fidelity. In Sec. IV, we consider quantum state 
teleportation by the Braunstein and Kimble protocol , 
and we identify the optimum operation given the input 
state distribution, and we obtain analytical results for the 
teleportation fidelity. In Sec. V, we consider a recently 
implemented quantum memory protocol 0, and we also 



here identify the optimum operation and fidelities. In 
Sec. VI, we consider teleportation of non-Gaussian states 
as exemplified by a distribution of Fock states. Section 
VII concludes the paper. 



II. GENERAL TRANSFORMATION OF AN 
UNKNOWN QUANTUM STATE 

Consider a family of states {|\I> m (A))}, parameterized 
by a stochastic variable A and used as input to a cer- 
tain quantum protocol according to the probability dis- 
tribution -P(A). For convenience, we assume the different 
states be obtained from one reference state |^o) by means 
of a family of unitary operators, \^> m (X)) = U\\^>q). 
The randomness of A is now handled by introducing an 
auxiliary (fictitious) physical system with no free evo- 
lution and with an initial mixed state, P (A)|A) (A|, 
where we assume that the quantum states |A) are or- 
thonormal. In our theoretical modelling, we let the 
auxiliary system interact in a Quantum-Non-Dcmolition 
(QND) manner with our physical system prepared in 
l^o), H = J2x \^}(M^H\, which we assume will lead, af- 
ter a suitable interaction time, to the following correlated 
state, 



p = £p(A)|A)<A|®|* in (A)) (¥ n (X)\. 



(1) 



The partial trace over the auxiliary A- 
degrees of freedom produces a density operator 
J2x p ( X W m ( X )) (* m (A)| describing the mixture of 
input states to the protocol with the appropriate proba- 
bilities, but we note that the fidelity of a protocol is not 
a measure of how well a mixed state is transformed into 
its image by the ideal protocol, but a measure of how 
well each member of the mixture transforms. It should 
also be noted that Eq.Q) identifies specific individual 
pure state components in the input, whereas the reduced 
density matrix £ A P(A)|* m (A)) (* m (A)| does not have 
a unique unravelling in terms of pure states. In |Q we 
have retained the variable A in the auxiliary system, 
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which enables us to study the fidelity at the pure state 
level, without selecting a specific pure state input to the 
protocol. 

A physical transformation of a quantum state must be 
completely positive and preserve the normalization of the 
density operator, and it can be written most generally in 
the Kraus form, p — > ^ s E s pEJ, where the E s operators 
can be any set of operators that fulfils J2 S -^1^ S = 1- 
Important examples include (i) unitary evolution, where 
there is only one unitary operator, E s = U, (ii) open sys- 
tem dynamics following a Lindblad form master equation 
with jump and no-jump operators E s , (iii) von Neumann 
measurements of a hermitian system observable with or- 
thogonal projections E s , and (iv) more general measure- 
ment positive operator valued measure (POVM) scenar- 
ios. 

Without loss of generality we assume that the desired 
protocol takes our quantum state to a final state on a 
similar Hilbert space (same dimensionality), so that each 
input state |^ m (A)) is ideally transformed by a unitary 
operation V into \^> out (X)) = VU\\^ ). We assume that 
physical interaction and measurements take place only 
on our quantum system of interest and on possibly added 
quantum systems.The variable A is not made subject to 
interaction or direct observation, and we hence write the 
state after application of the protocol, tracing over pos- 
sible further unobserved quantum degrees of freedom as 



£P(A)|A)(A|® \Y,E s \W n {X)){* in (ml) ■ (2) 



The sum over E s terms may represent the result of 
measurements on the system and it may also incor- 
porate a unitary feedback [/f° cdback on the system, 
conditioned on the outcome s of the measurement. 
(E s — ► iff eedhack E s a i so fulfil the required property 

^ g (f7feedback £ . s )t[/·feedback jE ; s ) = j Q f the Kraus form ) 

Equation J2J) is very illustrative. The sum over the dif- 
ferent Kraus operators corresponds to averaging over the 
outcomes of measurements and potential feedback on the 
system. It shows how each input component transforms 
into a mixed state p out (X) = J2 S E s \^ in (X)){^f in (X)\El. 
If we accept the output for all such measurement results, 
we must carry out this average, and we use the state in 
Eq. to compute the fidelity of the protocol. To check if 
p OMÍ (A) equals the desired state \^ out (X)) = VU\\^ in (X)) 
we apply the unitary J2x \^)(M ® (VUx)^ 1 on J2Jl, and 
check if the quantum system is now in the initial refer- 
ence state l^o)- Rather than verifying that the output 
equals the ideally transformed input state, we check if 
the inverse of the tranform on the output agrees with 
the fixed reference input state. Ideally, this agreement 
should be obtained for all A-components of the system 
and we thus perform the partial trace over the A-degrees 
of freedom and compare the final state system density 



matrix with the pure state \^o): 



F = J2PW(yv\UÍVÏp° ut (\)VU x \y Q )- (3) 



In this scheme we compute the average value of the over- 
lap without having to specify which initial state is ap- 
plied. 

We observe that this result can also be written 

F = ^P(A)(* ouí (A)|p ouí (A)|* ouí (A)), (4) 



where |* ouí (A)) = V\^ m (X)} is the desired output state, 
and despite our concerns in the Introduction about a con- 
sistent treatment of unknown input states, the fidelity is 
simply the fidelities obtained for each input state aver- 
aged over the input state distribution. A measurement 
part of the protocol may yield some information about 
the A- variable and hence change the probability distribu- 
tion P(A), but when we average over the outcomes s, we 
return to the original distribution. The fact that a feed- 
back may be applied to the system conditioned on the 
measurement affccts only the fidelity through the form of 
the E s operators. We note that Eq.(£j} could give the er- 
roneous impression that the P(X) distribution only enters 
via the explicit weighted sum. This is not the case; as we 
shall see below, to yield the highest possible fidelity the 
optimum feedback, i.e., the operators (/f ccdback should 
be chosen in a manner that depends on the distribution 

If the tranformation can be applied with a non-unit 
success probability, i.e., if the output state is only ac- 
cepted conditioned on a specific outcome or set of out- 
comes {s'} of the measurement on the system, we must 
go back to the joint state restrict the sum to only 
these vàlues and renormalize the state. (The trace of 
the un-normalized state is precisely the probability of 
acceptance). The resulting state is now a weighted 
sum of density operator terms with non-unit trace q\ = 
Tr (E s < E sl \^ m (X))(^ m {X)\El, ). We can renormalize the 
density operators with 1/qx and multiply the same q\ 
factors on P{X) which represents then the updated prob- 
ability distribution of the input states conditioned on the 
measurement result. Our fidelity calculation proceeds as 
above with the inverse transformations and the final com- 
parison with the initial reference state, and in this case, 
the fidelity is again given by the state-to-state transfor- 
mation fidelities but now weighted by both their initial 
state probability and their individual probabilities q\ for 
the acceptable measurement outcome. 



III. APPLICATION TO CONTINUOUS 
VARIABLE SYSTEMS 

Quantum information protocols with continuous vari- 
able systems have been the focus of intense research since 
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it was suggested [l| and demonstrated |3j that exist- 
ing squeezed light sources, beam splitters and photode- 
tectors suffice to enable quantum state teleportation of 
light. The collective atòmic population of different in- 
ternal states in a macroscopic gas sample also providcs 
cffcctively continuous degrees of freedom, and efficicnt 
atòmic entanglement protocols that make use of classi- 
cal light sources and photodetection only were proposed 
Q| and demonstrated . The work on entangled atòmic 
gasses was followed by theoretical and experimental work 
on quantum state transfer between light and matter (a 
quantum memory for light 0, IÜ ) , and ideas for atòmic 
state teleportation jj] are currently being pursued. 

These continuous variable systems can be described in 
terms of canonically conjugate harmònic oscillator vari- 
ables x and p, and states can be described in terms of 
Wigner phase space distribution functions in place of 
the general density matrix notation of the previous sec- 
tion. We consider the case where the ensemble of in- 
put states is obtained by displacements of the reference 
state (the vacuum state in Secs. IV and V) by argu- 
ments x c i and Pci according to a probability distribu- 
tion P(X — (x c i,pci)). Such a displacement of a Wigner 
function simply amounts to a translation of its argument 
W(x,p) — > W(x — x c \,p — Pel), but as in Sec. II we shall 
introduce an auxiliary set of QND variables in the mod- 
clling of the input state ensemble. We thus treat the real 
arguments x c \ and p c \ as two independent variables, e.g., 
momenta for free partides, or simply as classical vari- 
ables in a quantum-classical hybrid Wigner function for 
the total system, which is consistent with Heisenberg's 
uncertainty relation for the quantum degrees of freedom, 
but has no such constraints on the classical degrees of 
freedom. If we take the zero amplitude coherent state 
with a Wigner function Wo(x,p) and displace it by the 
classical arguments x c i,p c i according to a classical prob- 
ability distribution P(x c \,p c i), the joint Wigner function 
of the quantum and classical variables become in analogy 
with (Q, 

W- m {x,p,x c \,p c \) = W (x - x ch p - p c i)P(x ch p c i). (5) 

Some quantum information protocols make use of ad- 
ditional quantum systems and we shall hence work with 
a multi- variable Wigner function for all the quantum sys- 
tems and classical variables involved in the protocol. In 
teleportation, for example, the communication channcl 
is described by a joint Wigner function of the entan- 
gled state of two quantum systems W on t(xi,pi, X2,P2)- 
The total Wigner function is thus a function of 8 vari- 
ables W(x,p, x c \,p c \, x\,p\,X2,P2)- The linear mode mix- 
ing transformations of the teleportation protocol amount 
to the application of linear transformations on the vari- 
ables within the original distribution function; measure- 
ments of a given quantum variable amounts to evaluating 
the function with the corresponding argument attaining 
the measured value (and integrating over the canonical 
conjugate variable which is accordingly completely un- 
determined), and finally a joint distribution of the out- 



put quantum state and the classical variables is obtained. 
The verification of the protocol consists in comparing the 
output state with the desired one (which for teleporta- 
tion is the same as the input state) and this is done by 
displacing the quantum system with the negative of the 
classical variables (inverse of Eq.JSJ)), integrating over the 
unknown classical variables and comparing the ensuing 
quantum state with the reference state Wo(x,p). 

In the most general case one has to deal with a multi- 
variablc function, and one has to carry out integrals of 
this function with respect to a number of the variables. 
We shall now turn to examples where the initial quan- 
tum states and the classical distribution function are all 
Gaussian. This situation is of practical relevance in the 
above mentioned experiments and it offers a significant 
simplification of the problem. Gaussian states are fully 
characterized by the mean vàlues and the covariance ma- 
trix for the variables, and quantum state overlap inte- 
grals are given explicitly by these quantities. The present 
approach to the fidelity problem, involving joint prob- 
abilities for classical and quantum variables, is related 
to recent applications of the general theory of parame- 
ter estimation and Kalman filtcring with quantum sys- 
tems 0, ll| , and in particular to our recent application of 
this theory to atòmic magnetometry 0, . The formal 
treatment of the fidelity issue is actually simpler than the 
magnetometry analysis, and we now present the details 
of such a calculation. 



IV. TELEPORTATION OF AN UNKOWN 
COHERENT STATE 

We treat the case of teleportation of a physical system 
3 by use of an entangled pair of systems, 1 and 2. It has 
been argued, that a general positive map can be viewed as 
"teleporting a state through a gate" \LÏ\. and hence this 
operation has both specific and more general interest. 

We note that for coherent states with an amplitude of a 
given absolute value but with a random choice of complex 
phases, Ide et al. studied how to optimize the telepor- 
tation fidelity by a proper choice of the strength of the 
feedback on the output fields ■ We will consider states 
with a Gaussian distribution of complex amplitudes, and 
we will also identify an optimum feedback strength. Fi- 
urasek has applied the covariance matrix formalism, sim- 
ilar in spirit to our work, but rather than optimizing he 
assumed a fixed value for the feedback strength, and then 
he turned to a study of the effect of further local Gaussian 
operations 

A. Covariance matrix method 

We shall be working with Gaussian states, and hence 
the state is fully characterized by the mean vàlues m, = 
(yi) and the covariance matrix 7y = 2Re(((yj — {yi)){yj — 
(yj)) of all quadrature variables j/í. For a more detailed 
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description of the covariance matrix formalism and its 
practical implementation of linear transformations and 
measurementprocesses, see, e.g., [Ï4j . 

Following [lj, we use the entanglement in the 12-system 
to teleport an unknown coherent state of system 3 drawn 
from an ensemble of states with a Gaussian distribution 
of the mean amplitude onto system 1 by performing dis- 
placements of system 1 conditioned on the output of joint 
measurements on systems 2 and 3. We introduce the 
auxiliary classical variables {x c \,p c \) with zero mean and 
variance given by v c = 2Var(x c i) = 2Var(p c i). The clas- 
sical variables represent an agent Victor who's role is to 
turn the vacuum input state of variables (xs,pa) into a 
coherent state by a displacement. Experimentally, one 
is interested in the case with sizable v c (truly unknown 
input states), but we shall obtain general expressions for 
arbitrary v c . Note that since v c describes classical vari- 
ables, it is not limited by the Heisenberg uncertainty re- 
lation, and v c = corresponds to the case, where the 
input is the vacuum state with certainty. The covariance 
matrix of system 3 and Victors classical variables prior 
to the displacement is j 3 y = diag(l, 1, v c , v c ), and the 
displacement leads to the transformation 73^ — » j' 3 v = 

Sd 73, v , where the matrix Sd describes the linear map- 
ping x 3 — > x 3 +x c \ , p 3 -* P3 +Pd , Xd — > Zci, Pd -»■ Pel- The 
mean vàlues are also formally transformed by this map- 
ping, but since the classical distribution and the quantum 
state are assumed to have vanishing mean vàlues initially, 
this is also the case after the action of Sd- 

After this preparation of a quantum input state cor- 
related with classical stochastic variables as described 
by 73 v , we obtain the complete 8x8 covariance ma- 
trix 7 = blockdiag(7i2, 7 3 y), with 712 the covariance 
matrix for the initially entangled 12 system. The 
commuting pair of non-local variables x^ 23 ' = (X2 — 

x 3 )/\/2, p+ 3 ^ = (p2 + pa)/V% is measured. It is use- 
ful to transform the system to the following eight vari- 
ables: (xi,í»i, x 1 ^ 3 ^ ,p+ 3 \ x^ 3 \p^ 3 \ a; c i,Pci) with the co- 
variance matrix 7 — ► 7' = TjT T where the block diago- 
nal matrix T = blockdiag(/ 2 , T^, I 2 ), h the 2 x 2 iden- 
tity matrix (system 1 and the classical displacements are 
not affected by this transformation), and T^ 23 ) the matrix 
effecting the change of basis from system 2 and 3 variables 
to the joint variables x^ 3 ^ , p± . A measurement of a sin- 
gle variable from a joint Gaussian distribution results in 
an updated Gaussian distribution for the remaining un- 
known variables. This update is readily accounted for in 
terms of the mean vàlues and the covariance matrix of 
the variables. First, we reorganize the variables in the 

j / (23) (23) (23) (23) n , 

order (xi,pi,x c \,p c \,x + ,p + ,x_ ,p_ ), so that the 
covariance matrix 7' is decomposed into 4x4 dimensional 
matrix blocks: 

where A is the covariance matrix for the unmeasured 
quantum and classical components, B is the covariance 



matrix for the measured variables, and C and C T de- 
scribe the correlations between the measured an unmea- 
sured variables. The effect of the measurement on a sub- 
system on the covariance matrix of the remaining, un- 
measured variables is given by the update formula 0,0 

A' = A-C{-kBtt)-C t , (7) 

where 7r = diag(0, 1, 1,0) with unity at the entrances of 
P+ 3 \ and (. . . )~ denotes the Moore-Penrose pseu- 
doinverse. 

The measurement outeome affects the mean vàlues: 

conditioned on a positive readout £23 in the measure- 
(23) 

ment of x_ our knowledge about x 3 and hence of x c \ is 
biased towards negative vàlues and X2 is biased towards 
positive vàlues (and hence x\ is biased towards positive 
vàlues). Precisely how much, is determined by the vari- 
ances initially ascribed to these variables, and we have 
the following formula [Ï4l IÏ5j for the vector of mean vàl- 
ues, m= ((xi,pi,x ch p c i) T ): 

m^C(irBn)-(;r)23,& 3 ,-) T , (8) 

where the dots replace unmeasured quantities, which do 
not need to be specified due to the zeros in the projec- 
tor 7r, and 7723 and Ç23 denote the random outeome of 
the measurements of p+ 3 ^ and x^ 3 ^ . In Sec IV B, we 
shall give the expressions for the change of mean vàlues 
conditioned on the random measurement outeome. It is 
optimal to apply a feedback that brings the mean value 
of the quantum variables in exact agreement with the 
mean value of the classical displacement terms. This is 
so, because the covariance matrix does not depend on 
the measurement outeome, hence at the end of the cal- 
culation we shall compare the vacuum Gaussian state 
with either a single Gaussian state with vanishing mean 
or with a distribution of Gaussian states with the same 
covariance matrix but with different displacements with 
respect to the desired state. 

The feedback, just described is part of the teleporta- 
tion protocol. As part of our verification or quality as- 
sessment of the protocol, we displace the final quantum 
state by the negative of the classical parameters and com- 
pare the outeome, known to have vanishing mean ampli- 
tudes, with the vacuum state. Correspondingly, we apply 
the inverse of the classical displacement S^ 1 on the quan- 
tum and classical variables x\,p±, x c \,p c i and obtain their 
resulting covariance matrix, V — S d ~ 1 A'(S~ 1 ) T . The 
2x2 block 7 out = V(l : 2, 1 : 2) describing the covari- 
ances of the quantum variables is the covariance matrix 
for the quantum system, when the unknown classical dis- 
placements are integrated out, and it should ideally be 
the identity matrix describing the vacuum state. 

The fidelity of the protocol is the overlap of the Wigner 
functions. For a single mode state the Wigner function 
is given by W — l/(7T\/det7) exp(— x T 7 _1 x) with x T = 
(x,p). In terms of Wigner functions, the average fidelity 
is defined as F = 2iv j^' oo dx jT 00 dpW ín (x,p)W ut{x,p). 
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The vacuum state is Gaussian with a covariance matrix 
equal to the identity I 2 , and the integrant is thus a Gaus- 
sian function oc exp(~x T 7rcsX) wit h Tres = (7out + hY 1 
so that the integral follows directly from the Standard 
expressions for Gaussian normalization integrals, 



F 



I det(7 rcs ) 
det(7out) ' 



(9) 



B. Results 

For simplicity we consider the symmetric case wherc 
the joint covariance matrix of the variables (xi,pi, X2,p 2 ) 
for systems 1 and 2 is given by 



712 



/ n k 
n -k 
k n 
V -k n 



(10) 



where n describes twice the variance of the variables of 
system 1 and 2, and where k describe the correlations 
between the systems. The collective variables (x\ ± £2) 
and (jpi have the variances (n± k), and the Heisen- 
berg uncertainty relation implies that n 2 — k 2 > 1. Re- 
alizations of such a bipartite entangled state include the 
atom-light setup || and the EPR-light source channel Q . 
The matrix operations just described are straightforward, 
and we readily obtain analytical results at all steps of 
the calculation. The measurement process yields ran- 
dom outcomes £231 ?723; an d inserting the initial covari- 
ance matrices described above and carrying out the ma- 
trix operations, we obtain the conditioned mean vàlues 



(V26s) and 



similar expression for p\,p c \ in terms of the measured 
quantity 7723. At this point in the teleportation protocol, 
the aim is to have a state with (x\) = (x c \), (pi) = (p c i), 
and this is obtained by applying a feed-back on system 
1, in form of a displacement for both the x\ and pi vari- 
ables: 



Ífc(V2& 3 ), and 



x\ — > x\ 



Pi -+Pi 



k + v c 
1 + n + v c 

k + v c 
1 + n + v c 



(V263), 



(11) 



In the limit of infinite v c the feedbacks (|llf> are \Í2 times 
the measured vàlues themselves, but for states chosen 
from a finite width distribution, we see that a non-trivial 
gain factor 
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k + v c 



1 + n + v c 



(12) 



should be applied in the feedback. 

The resulting explicit expression for the fidelity reads 



F 



2{l + n + v c ) 



(1 + 2n + n 2 - k 2 + 2v c (l + n - k))' 



(13) 



If the input state is the vacuum state with certainty, v c = 
0, according to Ijl2|l . the optimum feedback gain is g = 
k/(l + n), and we note that for n 2 — k 2 = 1, which 
characterizes a pure two-mode squeezed state, system 1 
is restored in the vacuum state with unit fidelity. 
In the opposite limit v c — > 00, F simplifies to the result 



F 



1 



(14) 



where A = n — k = Var(xi — x 2 ) = Var(pi + po) is also 
known as the EPR variance of systems 1 and 2 [ïg . The 
fidelity approaches unity when this variance approaches 
zero corresponding to a maximally entangled channel. 
Finally, if the quantum channel is in the vacuum state 
with n = 1 and k = 0, the optimum gain is g — v c /(l + 
ri + v c ), and our general relation (|T3)) reduces to F = 
(2 + v c )/{2 + 2v c ) = (1 + A)/(2 + A) where A = 2/v c , 
which is exactly the best result that can be obtained 
with a classical strategy • 



V. QUANTUM STATE STORAGE 

Let us now turn to another example: an atòmic quan- 
tum memory, as demonstrated in a recent experiment 0] • 
In this protocol the aim is to store the quantum state of 
a light pulse in the collective spin degrees of freedom of a 
spin polarized atòmic sample. The transverse quantum 
degrees of freedom of the collective spin can be effectively 
described by canonical conjugate variables (their commu- 
tator, the polarized spin component, can be treated as a 
constant). In the protocol investigated in @], the optical 
Faraday rotation provides the light-atom interaction, de- 
scribed by a bilinear interaction Hamiltonian oc plPa- In 
the Heisenberg interaction picture this causes a change 
in the conjugate variables 



xa ■ 

XL 



x A - 

* XL 



kPl, 

- Kp A , 



(15) 



where k is the dimensionless integrated interaction 
strength 2] . The interaction thus encodes the field vari- 
able pl onto the atòmic xa, and by subsequently de- 
tecting the xl component of the field and displacing pa 
according to the measurement result, also this field com- 
ponent is read onto the àtoms. 

A theoretical analysis of the fidelity of this approach, 
applied to an unknown coherent state of light taken from 
a Gaussian distribution of coherent state amplitudes fol- 
lows the above discussion of teleportation. We introduce 
classical variables a; c i,p c i with variance parameter v c and 
zero mean and quantum variables for the àtoms and light 
in zero mean field coherent initial states, so that the 
Wigner function is a function of six variables, and the co- 
variance matrix is 6 x 6. We apply the linear transforma- 
tion between the field variables and classical variables to 
initialize the ensemble, and we apply the time evolution 
due to the atom-light interaction l|15|) . These operations 
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cause a mathematical transformation of the covariance 
matrix, and all mean vàlues still vanish. The detection 
of the field component Xl leads to an output value Ç, and 
conditioned on this output, we obtain the mean vàlues, 
(j) A ) = k/(1 + k 2 + v c )Ç and (x c i) = v c /(l + k 2 + v c )Ç. 
We wish to encode — xl in pa |2|, and shall hence apply 
a feedback on the atòmic pa variable p A — > Pa — ff£ with 
the non-trivial (optimal) gain factor 

9= (n + v c )/(l + k 2 +v c ) (16) 

As for teleportation, in the v c — > oo limit the optimum 



This general expression for the storage fidelity has sev- 
eral interesting limits. First, we observe that for a com- 
pletely unknown initial state with v c — > oo, the fidelity 
vanishe s un less K = 1, in which case one gets the value 
F = y/2/3 ~ 0.8165, reported in the literature 0- In 
the opposite limit of a known vacuum input, the choice 
k = 1 yields F = 2^2/3 ~ 0.9428 for the storage fidelity. 
The optimum strategy for finite v c , however, is to adjust 
the value of n so as to maximize the fidelity, and this 
leads to unit storage fidelity for n = 0, v c = 0. 

In the storage protocol, the field variables are mapped 
onto the atòmic ones, but part of the initial atòmic noise 
in the xa variable remains in the atòmic system whereas 
the feedback manages to cancel the pa component ex- 
actly. It has therefore been suggested to use an initially 



This expression yields the storage fidelity, and we again 
see that if v c — > oo, F vanishes unless k = 1. 



For infinite, or very large, v c and n = 1, the fidelity 
can be expanded in the squeezing parameter r <C v c , 
F — y/2r/(2r + 1) which yields the already mentioned 

y/2/3 for r — 1 and which approaches unity for large r. 



feedback gain is unity, but for finite width distributions 
it depends explicitly on the variance v c . The state stored 
is now guaranteed to have the same mean amplitudes as 
the classical variables, and to check if we managed to 
store pl in xa and — xl in pa, we follow the procedure 
from above and displace xa by ~p c \ and pa by x c \, and 
compare the resulting Gaussian state covariance matrix 
with the vacuum state as in Eq. (JÜJ. 

The resulting fidelity is a function of the variance of 
the classical variables and the coupling strength k: 



(17) 

I 

squeezed atòmic state. This is readily analyzed in our 
description. We simply take the vàlues (1/r, r) with r a 
squeezing parameter larger than unity for the initial di- 
agonal elements of the atòmic covariance matrix in the 
(xa , Pa) basis and go through all of the above steps again. 
In this case, we find the optimum feedback gain factor 

g = (nr + v c )/(l + n 2 r + v c ) (18) 

The state stored is again guaranteed to have the same 
mean amplitudes as the classical variables, and the fi- 
delity of the memory storage is a function of the initial 
atòmic squeezing, the variance of the classical variables 
and the coupling strength n: 



(19) 

I 

VI. TRANSFORMATION OF NON-GAUSSIAN 
STATES 



For quantum computing it is necessary to be able to 
store logical states and 1 and their superposition states, 
and one may not restrict an analysis to Gaussian states 
only, but still the above mentioned protocols may be use- 
ful. A Gaussian entangled state may be used to teleport 
also non-Gaussian states and a quantum memory 
protocol which transforms Gaussian states into Gaussian 
states also applies to qubit states, encoded in a two- 
dimensional subspace of the continuous variable Hilbert 



F 



(1 + V C K 2 — 2v c K + V c + K 2 + 2v c + 1)(1 + V C K 2 — 2v c K + V c + K 2 + 1) 



J 




r(l + K 2 r + v c ) 

(r + rv c K 2 — 2rv c K + rv c + K 2 r + 2v c + l)(r + rv c n 2 — 2rv c K + rv c + n 2 r + 1) 
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space @. 

A wide class of non-Gaussian states can be obtained 
by application of non-Gaussian operations on Gaussian 
states, e.g., photon counting on squeezed states [ï^ |. 
The corresponding non-Gaussian Wigner functions can, 
in turn, be expressed in terms of simple mathematical 
operations on Gaussian functions. This implies that re- 
sults explicitly derived for Gaussian states can be used 
as generating functions for quantities of relevance also for 
non-Gaussian states. Closer to the spirit of the present 
paper, we shall, however, give a few examples, where we 
explicitly apply the method of Sec. III, i.e., we assume 
that the Wigner functions are known for the input states 
to the protocol, and we carry transformations on the joint 
Wigner function of the entire physical system. 

We shall focus on teleportation by use of a two-mode 
squeezed state, and we shall apply the same protocol as 
above, but the input states will be number states and 
randomly displaced number states. 

The entanglement channcl of modes 1 and 2 is de- 
scribed by a covariance matrix 712, and hence by the 
two-mode Wigner function, 

W 12 {x 1 ,p 1 ,x 2 ,P2) = 7= exp(-x T 7r 2 1 x) ; (20) 

7r v /det(7i 2 ) 

where x T = Pii x 2iPi)- We assume the covariance 
matrix in Eq.(10), with the allowed vàlues of the param- 
eters n and k, listed in Sec. IV. B. 

The N — 1 Fock state input Wigner functions is given 

by 

W N ^{x 3 ,p 3 ) = -{2x\ + 2 PÍ ~ l)exp(-^ -pl). (21) 

TT 

The Wigner functions for the Fock states are all prod- 
uets of a polynomium in the arguments and a Gaussian. 
The beam splitter operations, the evaluation at the ar- 
guments measured, and the integration over conjugate 
and unmeasured variables, specified for the teleportation 
protocol in Sec. III, all preserve this mathematical form 
of the Wigner function, and it is hence possible to obtain 
analytically the outeome of the protocol and its fidelity. 
We shall now summarize the results of this analysis. 

We have evaluated the output state for different de- 
grees of entanglement of the teleportation channel. For 
N = Owe reproduce the results of Sec. IV. B with v c = 0, 
and for higher N we compare our results with Ref.|l8|. 
Our parameters n and k are equivalent to c and s in 
Ref.[Ï8j], and in the expression for the fidelity, our pa- 
rameter A = n — k is equivalent to t/2 in the notation 
of 0. (The expression t = 2/(c + s) in only applies 
for the pure state case, and should in the general case be 
replaced by 2(c — s) for the ensuing results to be correct). 
With these modifications, we reproduce the expression 
for the fidelities in |ïg|. and in particular the result 

1 + A 2 

Fn = i = otw (22) 



The analysis in 0] assumes a unit feedback. Applying 
instead a variable feedback gain g as in the previous 
sections, we are able to optimize the teleportation 
protocols also for non-Gaussian states: For the N = 1 
Fock state, we find that for strong entanglement (small 
A) unit gain is favored, but for weaker entanglement, it 
is advantageous to reduce the gain factor continuously 
to the value g = 1/^/2 for n = 1 and k = 0. This is 
summarized by the numbers 

{( n : ^iV=l,g=l) ^iV=l,aopt)} = 

{(8, 0.8364, 0.8524), (4, 0.7098, 0.7346), 

(2, 0.5258, 0.5602), (1,1/4 = 0.25,8/27 = 0.2963)} 

for pure state channels with n 2 = k 2 + 1. 

Teleportation of a known quantum state can in princi- 
ple be replaced by a local production of the given state 
with much higher fidelity. Let us therefore proceed with 
teleportation of unknown states, and let us begin with 
the teleportation of a Fock state, taken from an expo- 
nential distribution of Fock states, pn oc exp(— N/N), 
with mean value (N) — N. This distribution can also be 
written p N = (1 - X)X N with A = exp(-l/]V). The aver- 
age fidelity, according to our Eq. (4) is the mean value of 
the fidelities weighted with the probability distribution 
for the input states. In [ï^, the Fock state teleportation 
fidelities are derived from a generating function, which 
is, apart from a factor, precisely this mean value. We 
therefore readily obtain the mean fidelity as function of 
the channcl EPR variance and the parameter À 



y/(l + A) 2 - 2A(1 + A 2 ) + A 2 (l - A) 2 ' 

(23) 

The exponential distribution of N- vàlues is equivalent to 
a Gaussian distribution of amplitudes, hence the ensem- 
ble of Fock states has the same density matrix as the en- 
semble of displaced vacuum states treated in Sec. IV, but 
since we are dealing with the state-to-state teleportation 
fidelity, the results are very different. In particular, we 
found in (14), that when the variance of the distribution 
of coherent input states diverges, the fidelity approaches 
F = 1/(1 + A), whereas (12. '{t vanishes for fixed A and 
N = (v c — l)/2 — > 00. It requires a very strongly entan- 
gled channel, A < 1/A^, to reliably handle the difference 
between highly excited Fock states. 

In contrast, we have also implemented the scenario in 
which the input to the teleportation channel is the N = 
1 Fock state displaced by an unknown amount, similar 
to the displacements of the vacuum (N — 0) state in 
Sec. IV. With unit gain, in this case, we find that the 
result Í2'2\) holds irrespectively of the variance v c of the 
distribution of displacements. 



VII. CONCLUSION 

In conclusion, we have presentcd a theory to determino 
the fidelity of a general quantum state transformation on 
an unknown quantum system. The result of this analysis 
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is that, as long as the protocol has been definitely deter- 
mined in terms of the actions on the system conditioncd 
on the measurement outcomes, one can compute the fi- 
delity as a simple weighted average of the state-to-state 
fidelities over the incoming set of states. 

We have introduced a formalism which incorporates 
the preparation of the input states, and showed that our 
use of a fictitious system in a mixed state which is corre- 
lated with the input state to the protocol may indeed be 
convenient for practical calculations. We demonstrated 
this last point in the case of Gaussian transformations 
of Gaussian states, where we showed that the covariance 
matrix formalism readily identifies the optimum perfor- 
mance and provides simple analytical results for the fi- 
delity of teleportation and quantum memòries. The op- 



timal use of non-trivial feed-back gain in these protocols 
was a particularly interesting result brought out clearly 
by the analysis. Finally we showed that more general 
states can also be handlcd by their appropriate Wigncr 
functions. 
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